Proceedings Title IAU Symposium 
Proceedings IAU Symposium No. 241, 2007 
A.C. Editor, B.D. Editor & C.E. Editor, eds. 



© 2007 International Astronomical Union 
DOI: 00.0000/X000000000000000X 



Chemical clocks for early-type galaxies in 

clusters 

Conrado Carretero 1 , Alexandre Vazdekis 1 and John E. Beckman 1,2 

1 Instituto de Astrofisica de Canarias. 38200 La Laguna, Tenerife, Spain. Email: cch@iac.es 
2 Consejo Superior de Investigaciones Cientfficas, Spain. 

Abstract. We present a detailed stellar population analysis of 27 massive elliptical galaxies 
within 4 very rich clusters at redshift z~0.2. We obtained accurate estimates of the mean 
luminosity- weighted ages and relative abundances of CN, Mg and Fe as functions of the galaxy 
velocity dispersion, a. Our results are compatible with a scenario in which the stellar popula- 
tions of massive elliptical galaxies, independently of their environment and mass, had formation 
timescales shorter than ~1 Gyr. This result implies that massive elliptical galaxies have evolved 
passively since, at least, as long ago as z~2. For a given galaxy mass the duration of star for- 
mation is shorter in those galaxies belonging to more dense environments. Finally, we show that 
the abundance ratios [CN/Fe] and [Mg/Fe] are the key "chemical clocks" to infer the star for- 
mation history timescales in ellipticals. In particular, [Mg/Fe] provides an upper limit for those 
formation timescales, while [CN/Fe] apperars to be the most suitable parameter to resolve them 
in elliptical galaxies with <r<300 km s _1 . 
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C~-» , Stellar Population analysis has recently shown increasing power in constraining galaxy 
formation scenarios. They offer a fossil record of the star formation and chemical evolution 
of galaxies, most clearly in elliptical galaxies, thought to have formed by the merging 
Qh' of discs, so stellar population studies provide very strong constraints on these galaxy 
formation scenarios. For example, it is particularly hard to reconcile the hierarchical 
models with the result that massive galaxies show significantly older mean luminosity 
weighted ages than their smaller counterparts (Kauffmann et al. 2003). 

A key point is to see how abundance ratios in early-type galaxies vary with the cluster 
, masses. In particular, overabundances of [Mg/Fe] compared with the solar ratio have 
J^j ' been found in massive elliptical galaxies (e.g. Vazdekis et al. 1997). In this work, we 
concentrate on the stellar population analysis of 27 elliptical galaxies distributed in 4 
C$ ' rich, intermediate-redshift (z~0.2) Abell clusters. 



2. Observations and Data 

We have studied a set of 27 elliptical galaxies belonging to 4 very rich Abell clusters: 
A115, A655, A963 and A2111. The selected galaxies are the brightest galaxies in their 
respective clusters (15.1<my<17.2) and have been already morphologically analyzed in 
the literature. The velocity dispersion covers the range 220 km s _1 <er<450 km s _1 . 

Multi-slit spectroscopy was carried out with the DOLORES spectrograph on the 3.5m 
TNG Telescope (La Palma). We designed MOS masks with 1.1" width slits and an 
instrumental resolution of 8 A (FWHM). 
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Fe4383 CN2 Mg2 



Figure 1. Examples of the method used to derive ages and metallicities. The plot shows the 
age index H/? versus the metallicity indicators Fe4383, CN2 and Mg2 for the grid of models of 
Vazdekis et al. (2006, in preparation). Horizontal dotted lines represent models of constant age 
whose ages are labelled in the plot. Vertical dashed, solid, dotted, dashed-dotted, dashed-dot- 
ted-dotted and long-dashed lines represent contours of constant [M/H] =+0.2, 0.0, -0.4, -0.7, 
-1.2 and -1.7 dex, respectively. Diamonds, open circles and squares indicate the values for a rep- 
resentative galaxy of A115, A655 and A963, respectively, with a common <j=310 km s _1 . Note 
that all synthetic model spectra were broadened to the instrumental resolution of the observed 
spectra and the velocity dispersion of these particular galaxies (<7 o b s =350 km s _1 ). 

3. Analysis 

To derive mean luminosity-weighted ages and metallicities, we compared selected ab- 
sorption line strengths with those predicted by the model of Vazdekis et al. (2006, 
in preparation). This model provides flux-calibrated spectra in the wavelength range 
AA 3500-7500 A, at a resolution of 2 A (FWHM) for single-burst stellar populations. 

Once the model spectra are transformed to yield the intrumental conditions of resolu- 
tion and dispersion of the observed spectrum, we measure pairs of indices in both sets of 
data (observed and synthetic). Figure Q] illustrates this method. 



4. Results 

Figure [2] presents the relations [CN/H]-<7, [Mg/H]-cr, [Fe/H]-cr and Log(age)-er for our 
full sample. We have also plotted on Figure the results obtained by Sanchez-Blazquez 
et al. (2006), who analyzed 98 early- type galaxies of the Virgo and Coma clusters. 

Figure [3] presents the relations of the abundance ratios [CN/Fe] and [Mg/Fe] with the 
velocity dispersion. We observe positive trends with a for both ratios. 



5. Discussion 

First, when inspecting Figures [5] and we observe that, for a given cr, the relative 
abundance [CN/H] increases with the cluster X-ray luminosity, whereas [Mg/H] and 
[Fe/H] do not. The same happens with the abundance ratios: while [CN/Fe] increases as 
a function of the environment density, the [Mg/Fe] ratio keeps constant and positive. 

If we consider the formation timescales for Mg and Fe we can estimate that, whatever 
the environment is, all stellar populations in ellipticals should be already formed in less 
than ~l Gyr, since this is the time needed by SNe la to pollute the ISM with iron-peak 
elements. But, if we consider species with less separated formation timescales, such as Fe 
and CN, we do find differences in the ratio [CN/Fe] as a function of the environment, for 
a given a. That suggests a different formation timescale for the stars in ellipticals related 
to the environment properties, in the sense that the more dense the environment is the 
shorter the formation timescales of the stars are (with an upper limit of ~1 Gyr). This 
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Figure 2. Relations between the relative abundances [CN/H], [Mg/H], [Fe/H], and mean lu- 
minosity-weighted ages and the velocity dispersion of the galaxies, a. Diamonds: A115 galaxies. 
Open circles: A655 galaxies. Open squares: A963 galaxies. Open triangles: A2111 galaxies. The 
solid line corresponds to the error-weighted linear fitting to our data. Dotted line and dashed 
line corresponds to the relations obtained for Virgo cluster galaxies and Coma cluster galaxies, 
respectively. It is noteworthy that the slope of the relations increases with the environment 
density for [CN/HJ-u, while it is almost constant in the case of [Mg/H] -it and [Fe/H]-<r. 

is in agreeement with a number of recent results (e.g. Carretero et al. 2004, Thomas et 
al. 2005, Bernardi et al. 2006). 

Second, it is noteworthy that the slope of the relation [CN/Fe]-cr increases with the 
cluster X-ray luminosity while the slope of the [Mg/Fe]-<7 is constant with the cluster 
mass. Because of that, we can only find environment-related [CN/Fe] differences in those 
galaxies with cr<300 km s _1 whereas, for more massive galaxies, the [CN/Fe] values are 
almost equal independently of the cluster properties. This means that the ratio [CN/Fe] in 
most massive galaxies is less sensitive to the environment than in intermediate- and low- 
mass galaxies, suggesting that star formation histories are less environment dependent 
as the galaxy mass increases. Therefore, [CN/Fe] appears to be an appropriate "chemical 
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Figure 3. Abundance ratios [CN/Fe] and [Mg/Fe] versus the velocity dispersion of the galaxies. 
Symbols are the same than in previous Figure. Dashed-dotted line in bottom panel corresponds 
to the [Mg/Fe]-cr relation for Virgo cluster galaxies, by Vazdekis et al. (2004). Note that the 
slope of the [CN/Fe]-<r relation increases with the density of the environment, but the [Mg/Fe]-(j 
remains almost constant. 

clock" for ellipticals with er<300 km s~ x but, when we move to very massive galaxies, the 
formation timescales differ so gently from one environment to another that this "clock" 
turns out to be insufficiently sensitive. In order to disentangle the different formation 
timescales between very massive galaxies in different environments, it is neccesary to 
find out a more accurate indicator. 

Our results suggest an overall picture in which massive elliptical galaxies in very rich 
clusters are old systems, with very short star formation histories, and which have been 
passively evolving since, at least, z^2. This is in agreement with a number recent results 
(e.g. Labe et al. 2005, van Dokkum et al. 2006). In comparison with less dense environ- 
ments, the stars in elliptical galaxies of rich clusters must have formed at slightly earlier 
epochs and on a slightly shorter timescales. 
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Abstract. Using Stokes flow between eccentric counter-rotating cylinders as a prototype for 

bounded nearly parallel lubrication flow, we investigate the effect of a slender recirculation region 

C^l ■ within the flow field on cross-stream heat or mass transport in the important limit of high Peclet 

5^ ' number Pe where the 'enhancement' over pure conduction heat transfer without recirculation is 

most pronounced. The steady enhancement is estimated with a matched asymptotic expansion 

to resolve the diffusive boundary layers at the separatrices which bound the recirculation region. 

The enhancement over pure conduction is shown to vary as e 1 ^ 2 at infinite Pe, where e 1//2 is the 

characteristic width of the recirculation region. The enhancement decays from this asymptote 

as Pe~ 1/2 . 
t-H ■ 

^ ' Keywords. Keywordl, keyword2, keyword3, etc. 
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' 1. Introduction 

O ' 

The use of integral equations to solve 'exterior' problems in linear acoustics, i.e. to 
solve the Helmholtz equation (V 2 + k 2 )(f> — outside a surface S given that <fi satis- 
fies certain boundary conditions on S, is very common. A good description is provided 



O . by Martin (1980) Integral equations have also been used to solve the two-dimensional 
Helmholtz equation that arises in water-wave problems where there is a constant depth 



variation. The problem of wave oscillations in arbitrarily shaped harbours using such tech- 



o 

CZ2 - 

In a recent paper Linton & Evans (1992) have shown how radiation and scattering 



niques has been examined (see for example Hwang & Tuck 1970; Lee 1971; Figer, Najarro, Gilmore, et al. (2002) ) 



problems for vertical circular cylinders placed on the centreline of a channel of finite 
• i-h . water depth can be solved efficiently using the multipole method devised originally by 
X i |Ursell (1950)1 This method was also used by|Callan, Linton fc Evans (1991)]to prove the 



existence of trapped modes in the vicinity of such a cylinder at a discrete wavenumber 
k < n/2d where 2d is the channel width. 

Many water- wave/body interaction problems in which the body is a vertical cylinder 
with constant cross-section can be simplified by factoring out the depth dependence. 
Thus if the boundary conditions are homogeneous we can write the velocity potential 
<f)(x,y,z,t) — Re{0(x, y) cosh k{z + h)e~ luJt }, where the (x,?/)-plane corresponds to the 
undisturbed free surface and z is measured vertically upwards with z = —h the bottom 
of the channel. 

Subsequently Callan et al. (1991) proved the existence of, and computed the wavenum- 
bers for, the circular cross-section case. It should be noted however that experimental 

f Present address: Fluid Mech Inc., 24 The Street, Lagos, Nigeria. 
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evidence for acoustic resonances in the case of the circular cylinder is given by Bearman 
fc Graham (1980, pp. 231-232). 



Koch (1983) provided a theory for determining the trapped-mode frequencies for the 



thin plate, based on a modification of the Wiener-Hopf technique. Further interesting 



results can be found in Williams (1964) and Dennis (1985) 



The use of channel Green's functions allows the far-field behaviour to be computed in 
an extremely simple manner, whilst the integral equation constructed in §|3] enables the 
trapped modes to be computed in §|4] and the scattering of an incident plane wave to be 
solved in §[5j Appendix lAl contains comparisons with experiments. The Galactic Center 
region proves to be an area very rich in WR stars. The Vllth Catalogue lists within 50 pc 
from the Galactic Center 15 WNL and 11 WCL stars, at near-IR wavelengths discovered 
by Krabbe, Genzel, Eckart, et al. (1995) in the Galactic Center Cluster. 



2. Green's functions 

2.1. Construction of equations 

We are concerned with problems for which the solution, 0, is either symmetric or anti- 
symmetric about the centreline of the waveguide, y = 0. The first step is the construction 
of a symmetric and an antisymmetric Green's function, G S (P, Q) and G a (P, Q)- Thus we 
require 

(V 2 + k 2 )G s = (V 2 + k 2 )G a = (2.1) 
in the fluid, where V is a gradient operator, 

V-v = 0, V 2 P = V • (v x w). 



In ELL 



G s ,G a ~ l/(27t)lnr as r =|p_Q|-»o, (2.2) 
dG s 3G a 



dy dy 
dG s 



on y = d, (2.3) 



on y = 0, (2.4) 



dy 

G a — on y = 0, (2.5) 

and we require G s and G a to behave like outgoing waves as \x\ — > oo. 
One way of constructing G s or G a is to replace (|2.1[) and (|2.2p by 

(V 2 + k 2 )G s = (V 2 + k 2 )G a = -5{x - 0S(y - r,) (2.6) 

and to assume initially that k has a positive imaginary part. 

2.2. Further developments 



Using results from Linton & Evans (1992) we see that this has the integral representation 

/ ^[c-^y-^ +c- k ~< ( - 2d - y - r i')} cos k{x-£,)tdt 7 0<y, rj < d, (2.7) 

27T Jo 

where 

' -i(l- $2)1/2, i<l 



7(*) = 



(t 2 - 1)1/ 2 , t > 1. 
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In order to satisfy l|2.4p we add to this the function 

1 f°° N cosh fc 7 (d - y) ,. , 
-2^/ ^ 7 sinhfc 7 / W -Qtd* 

which satisfies (|2.ip . (|2.3p and obtain 

B(t) = 2e- fc T d coshfc 7 (d-?7). (2.8) 

Thus the function 



G 



1 f°° e~ kjd 

-\i{H {kr) + HoQen)) - / . , . , coshfc 7 (d - y) coshfc 7 (d - rf) (2.9) 

7T Jq 7 smn fc 7 d 

satisfies (|2.ip - (j2.4p . By writing this function as a single integral which is even in 7 , it 
follows that G is real. Similar ideas have been developed in a variety of ways (Keller 
1977; Rogallo 1981; van Wijngaarden 1968). 



3. The trapped-mode problem 

The unit normal from D to dD is n q = {-y'{9),x'{9))/w(9). Now G a = \Y {kr) + G a 

where r = {[x(9) — x(ip)] 2 + [%){&) ~ y(' l P)} 2 } 1 ^' 2 an d G a is regular as kr — > 0. In order to 
evaluate dG a (9,9)/dn q we note that 

as kr — > 0. Expanding x{ip) and y(tp) about the point ip = 9 then shows that 

' -[p'{9)p"{9) -p 2 {6) -2p' 2 {9)] as kr -> 0. (3.1) 



4nw 3 (9) 

3.1. Computation 

For computational purposes we discretize (13. ip by dividing the interval (0, n) into M 
segments. Thus we write 

M r) 

= ^^Wj)-fa G *W>A)M6j), 0<^<n, (3.2) 



where 9j — (j — ^)n/M. Collocating at ip = 9i and writing fa = <f>{9i) etc. gives 

M 

= Tl ^ h K >* ' i=h---,M, (3.3) 

3=1 

where 



j dG a {9 i ,9 :j )/dn q , i ^ j 

K * j = \ dGaWi, Oi)/dn q + Wl - Pf - 2p?]/4™,?, i = j. ' ' 4 * 

For a trapped mode, therefore, we require the determinant of the M x M matrix whose 
elements are 

2n 

Sij ~ Jl K "'"' r 

to be zero. Table Q] shows a comparison of results obtained from this method using two 
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Table 1. Values of kd at which trapped modes occur when p{6) — a 



n /W 
iljUL 


M — A 
IVl — *± 


IVl — o 


Cctllciii ct oil. 


0.1 


1.56905 


1.56 


1.56904 


0.3 


1.50484 


1.504 


1.50484 


0.55 


1.39128 


1.391 


1.39131 


0.7 


1.32281 


10.322 


1.32288 


0.913 


1.34479 


100.351 


1.35185 




WN subtype WC subtype 

Figure 1. Trapped-mode wavenumbers, kd, plotted against a/d for three ellipses: 
, b/a = 0.5; , b/a — 1; , b/a — 1.5. 

different truncation parameters with accurate values obtained using the method of Callan 
et al. (1991). 

An example of the results that are obtained from our method is given in figure [TJ 
Figure 0( a, b) shows shaded contour plots of cf> for these modes, normalized so that the 
maximum value of 4> on the body is 1. Symmetric (figure^) modes are shown, while 
the antisymmetric ones appear in figure^ b). 

3.2. Basic properties 

Let 

pi = lim p(x,C), p u = lim p{x,() (3.5a, 6) 

be the fluid densities immediately below and above the cat's-eyes. Finally let po and Nq 
be the constant values of the density and the vorticity inside the cat's-eyes, so that 

(p(x,(),(j> a (x,0) = (po,N Q ) for Z,(x) <t<Z u (x). (3.6) 

The Reynolds number Re is defined by u T H/v (y is the kinematic viscosity), the length 
given in wall units is denoted by ( )+, and the Prandtl number Pr is set equal to 0.7. In 
(j2~Tj) and (|2~2"| . Tij and r| are 

Uj = (uiUj - UiUj) + {uiU S } GS + uf GS u 3 ) + uf GS uf GS , (3.7a) 

7-/ = (^1 - U~6) + (UjOSGS + u SGSty + u SGS e SGS. (3.75) 

3.2.1. Calculation of the terms 

The first terms in the right-hand side of (|3.5b ) and (|3.5& ) are the Leonard terms 
explicitly calculated by applying the Gaussian filter in the x- and z-directions in the 
Fourier space. 

The interface boundary conditions given by (|2.1[) and (|2.2[) . which relate the displace- 
ment and stress state of the wall at the mean interface to the disturbance quantities of the 
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Figure 2. Shaded contour plots of the potential (f> for the two trapped modes that exist for an 
ellipse with a/d = 1.5, b/d — 0.75. (a) Symmetric about x = 0, kd = 0.96; (6) antisymmetric 
about x — 0, kd = 1.398. 



flow, can also be reformulated in terms of the transformed quantities. The transformed 
boundary conditions are summarized below in a matrix form that is convenient for the 
subsequent development of the theory: 



Qc = 



V 

2 w 



-(c^w)- 1 






(in') X>. (Si^+cW 1 ) 



iiVV+^O 







WRs)- 1 







i 









(3.8) 



S l is termed the displacement-stress vector and Qc the flow-wall coupling matrix. Sub- 



script w in (|3.8p denotes evaluation of the terms at the mean interface. It is noted that 
V", = for the Blasius mean flow. 



3.2.2. Wave propagation in anisotropic compliant layers 

From (|2.1[) . the fundamental wave solutions to (|3.1[) and (|3.2[) for a uniformly thick 
homogeneous layer in the transformed variables has the form of 

if exp[i(a*x* - ut)], 



ri 



(3.9) 



where if = bexp^a:^). For a non-trivial wave, the substitution of (1.7) into (1.3) and 
(1.4) yields the following determinantal equation: 



Det [pui 2 S ps 



(I 1 1 



0. 



(3.10) 



where the wavenumbers (k{, k\) = (a*, 0, 7) and S ps is the Kronecker delta. 



4. Torus translating along an axis of symmetry 

Consider a torus with axes a, b (see figure [2]), moving along the z-axis. Symmetry 
considerations imply the following form for the stress function, given in body coordinates: 

/(M) = (ffty)cosM^)siii0,/(VO). (4- 1 ) 
Because of symmetry, one can integrate analytically in the ^-direction obtaining a pair 
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of equations for the coefficients /, <? in (|4. 1 [) 

S.l/2 



7t[2(a + focosVi)] 3 / 2 i, 
3 



ln (sin^i - smtp)(a + 6cos^) 1/2 

[1 -cos(V>i - V)](2 + a)V2 ' ( ' > 



r 2n 



1/2 

= a zr / M3/2 / rt^ J {/W)[(cob^-6A)5 + AP] 

n[2(a + bcosipi)\ 6 l z J \ 2 + a J { 

sin^i -simp [ / (sin^ -sin^) 2 2 \ 
x i r; + ffW 2 + a - — — — 676 

i-cos(v>i-v) IA 1 - cos(^ - Vi) / 

+ (V cos ^17- ^aj F(±n,5) - (2 + a) cos foE^n, 6) \dip, 

a ,s 6 2 [1 -cos(^> - Vi)] o or 1 1, \ l-cos(V>-Vi) /, ,s 

a = a(V>,V'i) = 7 — ! r 7T7 — —r tt5 P ~ PW,i>V = tt ; — ■ ( 4 - 4 ) 

(a + cos -0) (a + b cos V>i J a + b cos V> 



(4.3) 



5. Conclusions 

We have shown how integral equations can be used to solve a particular class of prob- 
lems concerning obstacles in waveguides, namely the Neumann problem for bodies sym- 
metric about the centreline of a channel, and two such problems were considered in 
detail. 



Appendix A. Boundary conditions 

It is convenient for numerical purposes to change the independent variable in (|4.ip to 
z = y/v^ 2 and to introduce the dependent variable H(z) — (/ — y)/v^ 2 . Equation (|4.ip 
then becomes 



(1 - (3){H + z)H" - (2 + H')H' = H'". 
Boundary conditions to (I4.3[) follow from (|4.2p and the definition of H 

eC v 



H(0) = 
H"{0) 



4 /2 (i - /?) 



; H'(0) = -l + e 2 ^C u + eC' u 
H'(oo) = 0. 



(Al) 



(A 2) 



Appendix B. 

A simple sufficient condition for the method of separation of variables to hold for the 
convection problem is derived. This criterion is then shown to be satisfied for the ansatz 
described by (3.27), thus justifying the approach used in §[3J The basic ingredient of our 
argument is contained in the following estimate for a Rayleigh-Ritz ratio: 

Lemma 1. Let f(z) be a trial function defined on [0, 1]. Let A\ denote the ground-state 
eigenvalue for —d 2 g/dz 2 = Ag, where g must satisfy ztdg/dz + ag = at z = 0, 1 for 
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some non-negative constant a. Then for any f that is not identically zero we have 



Before proving it, we note that the first inequality is the standard variational charac- 
terization for the eigenvalue A\. 

Corollary 1. Any non-zero trial function f which satisfies the boundary condition 
/(0) = /(l) = always satisfies 
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Discussion 

Massey: Im wondering if you have considered the expected intrinsic dispersion in ab- 
solute magnitude of WRs - if you consider the (large) mass range that becomes an 
early WN or late WC according to the evolutionary models, wouldnt you expect a large 
dispersion in M„? 

VAN DER Hucht: Indeed, we will be always left with some intrinsic scatter in M„ due to 
mass differences within the same spectral subtype. But in my opinion, the current large 
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dispersion is for a large fraction due to incertaintics of the adopted distances of open 
clusters and OB associations. 

Walborn: I think that the scatter in WNL absolute magnitudes is dominated by intrinsic 
spread rather than errors. In the LMC, one finds a range of 5 to nearly 8. This in 
turn likely reflects different formation channels: mass-transfer binaries, post-RSG, and 
extremely massive stars in giant H n regions. 

van der Hucht: As said above, there is likely to be intrinsic scatter. But, I wonder 
whether a scatter of 3 magnitudes perhaps reflects undetected multiplicity. 

Maiz-Apellaniz: I could not agree more with your comment on the need for an updated 
catalogue of O-type stars (as a follow up of that of Garmany et al. 1982). We are currently 
working on precisely that (see our poster, these Proceedings) and we will soon make it 
available. 

van der Hucht: Wonderful. 

Koenigsberger: Is the ratio WR/O-stars in clusters similar or different from this ratio 
for the field stars? 

VAN DER Hucht: I think it is different because incompleteness among field stars is even 
larger than that among cluster stars. But perhaps it should also be different because WR 
stars are older and could have drifted away from clusters, more than O-type stars. 

GlES: How many of the WR stars in your catalogue might be low mass objects? 

Walborn: Comment: PN central stars in the WR sample would be only [WC]. 

VAN DER Hucht: Among the WR stars in our Vllth Catalogue we doubt only one: 
WR109 (V617 Sgr), which is a peculiar object (not even a [WR] central star of a PN). 
All other stars in our cataloge are true massive Population I WR stars, and properly 
classified as such. We have not listed known Population II [WC] objects, as we did 
separately in our Vlth Catalogue (van der Hucht et al. 1981). [WN] objects are not 
known to exist, see the comment by Nolan. 

Zinnecker: Are all Galactic WR stars in open clusters and OB association or are there 
many WR stars in the field? 

van der Hucht: See the Vllth WR catalogue (van der Hucht 2001): of the listed 227 
Galactic WR stars, only 53 are in open culsters and OB associations, or believed to be. 
The other 184 are supposedly field stars. 



